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ABSTRACT 

This  report  examines  the  feasibility  of  predicting  rates  of  crack  growth  in 
rndld  steel  based  on  an  evaluation  of  the  crack  tip  dissipated  energy.  A  numer¬ 
ical  scheme  is  developed  to  provide  a  basis  for  the  computation  of  dissipated 
power  from  measurements  of  temperature  acquired  using  remote  infrared  de¬ 
tection.  An  experimental  study  is  then  described  where  measurements  of  crack- 
growth  rate  in  mild-steel  coupons  exposed  to  constant  amplitude  cyclic  loading 
are  related  to  calculations  of  power  dissipation.  A  simple  linear  relationship 
is  indicated,  with  testing  also  providing  some  evidence  of  geometry  indepen¬ 
dence,  a  key  finding  with  regard  to  the  possible  use  of  dissipated  energy  as  a 
fatigue  characterisation  parameter.  However,  further  testing  involving  a  more 
complex  loading  sequence  led  to  a  noticeable  change  in  relationship  indicat¬ 
ing  that  the  heat  dissipation  rate  is  unlikely  to  offer  a  simple  basis  for  the 
prediction  of  crack  growth  rate  under  all  loading  conditions. 
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Thermographic  Assessment  of  the  Heat  Dissipation 
Associated  with  Fatigue  Crack  Growth  in  Mild  Steel 


EXECUTIVE  SUMMARY 

The  continued  assurance  of  aircraft  structural  integrity  relies  in  part  on  the  accuracy  of 
predictions  of  the  rate  of  growth  of  cracks  in  structurally-critical  members  of  an  airframe. 
When  based  on  the  stress  intensity  factor  such  predictions  occasionally  fail,  often  as  a 
result  of  the  inability  of  the  analysis  to  accoimt  for  plastic  behaviour  at  the  crack  tip. 
In  response  to  this  shortfall,  alternative  analysis  techniques  have  been  developed  which 
better  accoimt  for  the  role  of  crack  tip  plasticity  in  the  crack  propagation  process. 

This  report  examines  the  feasibility  of  predicting  rates  of  crack  growth  in  mild  steel 
based  on  an  evaluation  of  the  crack  tip  dissipated  energy.  A  numerical  scheme  is  de¬ 
veloped  to  provide  a  basis  for  the  computation  of  dissipated  power  from  measurements 
of  temperature  acquired  using  remote  infrared  detection.  An  experimental  study  is  then 
described  where  measurements  of  crack-growth  rate  in  mild-steel  coupons  exposed  to  con¬ 
stant  amplitude  cyclic  loading  are  related  to  calculations  of  power  dissipation.  A  simple 
linear  relationship  is  indicated,  with  testing  also  providing  some  evidence  of  geometry  in¬ 
dependence,  a  key  finding  with  regard  to  the  possible  use  of  dissipated  energy  as  a  fatigue 
characterisation  parameter.  However,  fiirther  testing  involving  a  more  complex  loading 
sequence  led  to  a  noticeable  change  in  relationship  indicating  that  the  heat  dissipation 
rate  is  urdikely  to  offer  a  simple  basis  for  the  prediction  of  crack  growth  rate  under  all 
loading  conditions. 
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1  Introduction 


The  inability  of  the  stress  intensity  factor  to  account  for  some  crack  growth  behaviours 
[1,2,3]  is  often  attributed  to  its  failure  to  acommodate  plastic  behaviour  at  the  crack  tip. 
This  shortfall  has  driven  the  development  of  alternative  crack  characterisation  parameters 
based  on  a  more  explicit  consideration  of  crack  tip  plasticity  ([4]  [5]).  Importantly,  it  has 
also  served  to  underscore  the  key  role  that  plastic  flow  plays  as  a  precursor  to  both  the 
initiation  and  growth  of  fatigue  cracks  in  ductile  materials. 

It  is  now  weU  imderstood  that  the  energy  expended  in  deforming  a  ductile  mater¬ 
ial  is  consumed  in  the  production  of  heat  and/or  stored  in  processes  of  crystalline  re¬ 
arrangement  involving  the  formation  and  motion  of  dislocations  and  vacancies,  and  if  crack 
growth  occurs,  an  increase  in  surface  energy.  Energy  storage  is  thought  to  be  strongly  re¬ 
lated  to  fatigue  damage;  the  most  intuitive  connection  being  that  between  crack  growth 
and  an  increase  in  surface  energy.  Experimental  studies  over  the  past  few  decades  have 
gone  to  some  length  in  supporting  this  notion  (see  [6,7,8]),  but  the  task  of  experimen¬ 
tal  verification  has  been  fraught  by  the  difficulty  in  measuring  the  relevant  energy  in  a 
non-intrusive  manner.  Two  techniques  have  traditionally  been  applied  to  energy  measure¬ 
ment:  annealing  calorimetry  and  the  single-step  method  [9].  The  former  entails  the  use  of 
a  calorimeter  to  measure  the  energy  liberated  by  recrystallisation,  where  the  stage  in  the 
annealing  process  at  which  the  release  occurs  can  assist  in  identifying  the  energy-release 
mechanism.  However,  in  imposing  a  microstructural  change,  the  method  is  not  suited 
to  applications  where  the  aim  is  to  track  the  evolution  of  stored  energy  as  a  function  of 
exposure  to  cyclic  loading.  Better  adapted  to  this  task  is  the  single-step  method,  where, 
based  on  energy  conservation  principles,  the  stored  energy  is  inferred  from  measurement 
of  the  applied  work  and  dissipated  heat,  viz.. 


Ail  =  AW-  AQ  (1) 

where  A [7  is  the  rate  of  change  of  stored  energy,  AW  is  the  rate  of  change  of  supplied 
work,  and  AQ  is  the  rate  of  change  of  dissipated  heat.  The  applied  work  is  easily  deduced 
from  an  appropriately  measured  force-strain  curve.  Evaluation  of  dissipated  energy  is  an 
inherently  more  difficult  exercise  as  it  relies  on  the  precise  measurement  of  temperature  at 
carefuUy  chosen  locations  on  the  structure.  Traditionally,  the  use  of  miniature  thermistors 
or  thermocouples  for  this  task  has  necessitated  the  adoption  of  tenuous  assumptions  in 
the  energy  calculations,  usually  comprising  the  over-simplification  of  thermal  boundary 
conditions  and,  in  extreme  cases,  can  limit  investigations  to  circumstances  where  thermal 
equilibrimn  is  established  (see  [10]). 

Recent  advances  in  infrared  detection  technology  have  greatly  eased  the  burden  of  ac¬ 
quiring  temperature  measurements  at  the  requisite  precision  for  useful  insitu  assessments 
of  dissipated  energy.  This  report  outlines  an  appropriate  methodology  that  exploits  the 
relatively  recent  innovation  of  focal-plane  array  infrared-detection  technology.  A  numeri¬ 
cal  scheme  is  developed  which  facilitates  the  prediction  of  energy  dissipation  rates  imder 
transient  conditions  based  on  radiometric  measurement  of  the  surface  temperature  dis¬ 
tribution.  Although  the  scheme  provides  a  framework  for  advancing  research  into  the 
finks  between  stored  energy  and  fatigue  damage,  the  experimental  work  described  in  this 
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report  will  focus  only  on  the  measurement  of  dissipated  energy,  with  the  aim  of  explor¬ 
ing  the  prospect  of  predicting  crack-growth  rates  based  solely  on  remote  thermographic 
assessment  of  energy  dissipation. 


2  A  Computational  Scheme  for  Estimating 
Energy  Dissipation  from  Temperature 
Measurements 


The  formulation  given  is  an  expansion  to  two-dimensions  of  the  one-dimensional  ap¬ 
proach  previously  used  [9]  in  a  thermographic  investigation  of  energy  storage  in  a  cylin¬ 
drical  specimen  exposed  to  sustained  cyclic  loading.  The  starting  point  here  is  the  heat- 
diffusion  equation  in  two-dimensional  rectangular  coordinates, 


dt  ^  dip  pC’ 


(2) 


where  9  is  the  temperature  change,  q  is  the  power  density,  a  is  the  thermal  diffusivity, 
p  is  the  density  and  C  is  the  specific  heat.  The  objective  is  to  find  an  estimate  for  the 
power,  q,  corresponding  to  some  measured  temperature  response,  a  problem  conveniently 
expressed  in  terms  of  the  minimisation  of  an  objective  function  X, 


X^  jjjT-efdxdy,  (3) 

where  T  is  the  measured  temperature,  9  is  the  temperature  response  associated  with  the 
imknown  heat  source,  and  A  is  the  domain  over  which  the  solution  is  sought.  Rewriting  the 
objective  function,  {Eqn  3),  to  include  the  constraint  that  the  xmdetermined  temperature 
response,  9,  satisfies  the  heat  diffusion  equation  {Eqn  2),  yields 


x  =  JJ  {T-9pdxdy-XF, 

(4) 

^  ae  ,0^9  dH.  q 

^^dx^'^dy^^  pC’ 

(5) 

where  A  is  the  Lagrange  Multipher  and  F  is  the  constraint  equation.  In  approaching 
the  optimisation  problem  numerically,  the  objective  function,  Eqn  4,  is  more  appropriately 
expressed  in  a  discretised  form,  i.e.. 


hi 

where  the  subscripts  i  and  j  refer  to  nodal  positions  in  the  rectangular  coordinates  x 
and  y  respectively.  The  contraint,  F,  is  discretised  using  a  Crank-Nicolson  scheme. 
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For  the  sake  of  simplicity  and  clarity  in  the  following  derivation,  mesh  spacings  in  the 
X  and  y  directions  are  assumed  equivalent  {Ay  =  As).  Accordingly,  Fij  is  written  as. 


0/\„2 


(7) 


Here,  nAt  is  the  time,  K  is  the  thermal  conductivity,  and  q  is  the  averaged  power 
density  over  the  time  increment  At,  viz.. 


I  rt+m 

%S  =  Xth 

For  convenience,  Eqn  2  is  expressed  in  the  following  abbreviated  form, 


Fid  =  ^^3  -  +  ^i-1,3  +  +  ^l3-l)  -  ^i,3  -  C'i.i, 


where, 


(8) 


(9) 


A  = 
B  = 


.2Aa:2  . 

2Arc2 


K  ’ 

on— 1  I  an—1 


,2A®2 


Substituting  Eqn  9  into  Eqn  6  then  yields. 


in— 1 


(10) 

(11) 

(12) 


AT-lM-l 
i=2  j=2 

DiflFerentiating  X  with  respect  to  the  unknowns  6,  X  and  q,  results  in  the  following  set 
of  normal  equations. 


-IM-l 

^  X- 


Aa. 


dX 

d9ij 


N-lM-l 

=  E  E  9i,j)Sij  AXij  +  Aj-j-i  j  +  Aj—i  j  +  Aij+i  +  Ajj_i 

i=2  j=2 


=  0, 


(15) 
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dX 


N-lM-l 

=  E  E 

i=2  j=2 


(16) 


It  is  possible  to  show  that  the  relevant  stationary  value  is,  as  required,  a  minimum. 
Additionally,  the  linearity  of  the  objective  function  ensures  uniqueness  of  the  solution. 
Consider  now  the  solution  of  the  normal  equations  (Eqns  14-16),  which  are  expressed  in 
the  matrix  form. 


AX  =  Y,  (17) 

where  A  is  a  matrix  of  coefficients,  and  X  and  Y  are  vectors  of  the  tmknowns  and 
forcing  functions  respectively.  The  heat  source  is  expressed  as  a  function  of  x  and  y  with 
K  terms,  namely. 


n,3 


{D\,D2,Dz,  ■  •■,Dk) 


/  \ 


(18) 


Here,  Di  are  the  unknown  regressor  parameters  and  Pi  are  the  regressor  terms,  which 
are  described  later.  The  vectors  X  and  Y  may  then  be  written  as. 


62.2 

62.3 

€2,2  —  22,1  —  ri^2 

C2,3  —  2u,3 

^2,2 

-^2,3 

’ 

Y=‘ 

Cn-i,m-x  -  2V,m-i  -  Tn-i,m 
2^2,222,2 

2^2,322,3 

' 

Di 

D2 

2^Ar-l,M-l2V-l,M-l 

0 

Dk 

0 

The  matrix  A  is  written  in  the  following  compound  form. 
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(20) 


where  Mi,  M2  and  M3  are  matrices  and  Hj(  and  Vk  vectors  comprising  the 
regressor  functional  terms  defined  in  Eqn  18. 

It  is  worthwhile  remarking  on  the  structure  of  the  component  matrices.  Mi  simply 
comprises  coefficients  of  the  Crank-Nicolson  problem  defined  earher  and  requires  no  further 
comment.  M2  includes  conditions  imposed  to  reconcile  Eqn  15  with  the  absence  of  6  on 
the  boimdaries,  i  =  1,N  and  j  =  1,  M.  This  problem  is  resolved  by  expressing  6  in  terms 
of  known  values  of  T  and  defined  values  of  A.  As  an  example,  consider  the  boundary  at 
*  =  1-  Noting  that,  like  9,  A  is  defined  within  the  interior  domain  i  =  2, 3, . . . ,  iV  —  1,  j  = 
2,3,...,M  —  1,  manipulation  of  Eqn  i 7 yields  the  expression. 


-2(21, j  -  dij)  +  A2,j  =  0,  (21) 

from  which, 

0x,j  =  Tij  -  (22) 

Repeating  this  for  the  other  boundaries  gives, 


Onj  = 

= 

2i,i  -  2. 

^i,M  = 

2i,M  —  1 

and  similarly  for  the  corner  points. 

^1,1  = 

2l,l  -  A2,2, 

^l,M  = 

T\,M  -  A2,M-1) 

Gn,i  = 

I — 1 

1 

1 

1— ( 

0N,M  = 

Tn,m  —  Ajv-1,M- 

(23) 


(24) 
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M3  is  a  modified  identity  matrix,  where  the  diagonal  values  represent  the  weightings 
applied  to  the  experimental  data  points.  Should  the  need  arise  to  exclude  certain  points 
from  the  analysis,  this  can  be  achieved  by  simply  assigning  a  value  of  0  to  the  appropriate 
diagonal  element. 

As  stated  earlier,  the  heat  source  fimctional  terms  appear  in  the  vectors  V  and  H, 
which  for  the  fcth  term  may  be  written  as, 


Vt  =  -B^Pk  =  -B\ 


‘‘P2.2  \ 


(25) 


and, 


Hk  =  =  {’‘P2,2,’^P2,3,  .  .  .  ,’^Pn-i,m-i)  ■  (26) 

Selection  of  an  appropriate  functional  form  for  the  terms  should  reflect  the  aims 
of  the  planned  experimental  work.  It  is  important  to  emphasise  that  the  objective  in 
this  study  is  not  to  recover  the  precise  spatial  distribution  of  the  energy  dissipation,  but 
merely  to  measure  the  total  crack-tip  energy  dissipation.  In  this  limited  context,  a  simple 
low-order  form  should  suffice.  A  biquadratic  expression  was  chosen,  viz.. 


q{x,y)  =  Dix"^  +  D2y^  +  D^xy  -t-  D^x  -|-  D^y  +  £>6,  (27) 

which,  in  the  form  of  Eqn  18,  is  written  as, 


%j  =  (Di,  D2,Dz,  Di,  D5,  De) 


xy 

X 


[IJ 


hi 


(28) 


Testing  of  the  scheme  initially  involved  synthetic  data  produced  from  finite-difference 
modeling  of  an  appropriate  two-dimensional  heat-dissipation  problem.  Estimates  of  the 
total  power  were  obtained  from  several  sets  of  synthetic  data  each  augmented  with  a 
different  amount  of  Gaussian  noise.  It  was  found  that  good  estimates  were  produced  even 
in  the  presence  of  imrealistically  high  levels  of  random  noise,  providing  verification  of  the 
precision  and  robustness  of  the  proposed  scheme.  As  expected,  the  scheme  was  far  more 
sensitive  to  systematic  errors  in  the  data.  These  tests  also  revealed  that  computations 
quickly  become  laborious  with  increased  mesh  size,  due  to  the  growth  of  the  coefficient 
matrix.  Fortunately,  robustness  was  formd  to  be  excellent  even  with  moderately  sized 
meshes.  Accordingly,  a  mesh  of  10  x  10  was  chosen  for  all  subsequent  experimental  work. 
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2.1  Experimental  Validation 

An  experimental  test  case  was  devised  to  facilitate  the  practical  verification  of  the 
outlined  scheme.  Figure  1  shows  the  experimental  set-up,  the  salient  aspect  of  which  is  the 
use  of  a  Peltier  ceU  applied  to  the  back  face  of  an  aluminium  sample  to  simulate  the  effect 
of  an  internal  heat  source.  Temperatme  response  data  were  acquired  through  inspection  of 
the  front-face  with  a  focal-plane  array  infrared  imager  (described  in  more  detail  shortly). 
Several  different  Peltier-cell  power  settings  were  considered.  In  each  case,  before  acquiring 
data,  a  resting  period  of  5  minutes  was  allowed  to  accommodate  the  re-establishment  of 
thermal  equilibrium.  This  was  necessary  to  facilitate  the  independent  computation  of  the 
effective  power  dissipation,  achieved  by  integration  of  the  normal  temperature  gradient 
along  a  closed  contour  surrounding  the  zone  of  heat  production  (see  Figure  2).  At  steady 
state,  it  can  be  shown  that. 


JJ^q{x,y)dA  =  ^^-K^ds,  (29) 

where  d/ dn  denotes  differentiation  in  the  direction  of  the  outward  normal  to  the  circuit 
path  and  ds  is  an  infinitesimal  segment  along  the  contour  C.  Estimates  derived  from 
the  optimisation  scheme  were  found  to  agree  well  with  those  derived  through  contour 
integration.  The  comparison  in  Figure  3  for  7  different  nominal  power  settings  shows  a 
discrepancy  that,  in  the  worst  case,  amounts  to  less  than  two  standard-deviations  of  the 
respective  optimisation  result. 


3  Experimental  Procedure 

Single-edge  notched  specimens  were  manufactured  from  a  low-carbon  steel,  and  used 
in  the  as-recieved  condition.  Two  configurations  were  considered,  shown  in  Figure  4  as 
t3rpe-A  and  type-B.  The  fatigue  testing  involved  sinuisoidal  loading  at  a  fixed  frequency 
of  12.5  Hz,  an  R-ratio  of  0.1  and  a  mid-section  stress  amplitude  of  122  MPa.  As  already 
stated,  temperature  measmements  were  acquired  with  an  infrared  imager.  Its  salient 
feature  is  a  mechanically  cooled  512x512  Platinum  SiHcide  photonic  detector  with  a  rated 
noise  equivalent  temperature  detectivity  of  150  mK  and  a  fixed  frame  rate  of  25  Hz. 

3.1  Analysis  of  the  Cyclic  Plastic  Zone  Size 

Initial  experimental  work  focused  on  mapping  the  cyclic  plastic  zone  produced  during 
crack  growth  at  the  stipulated  loading  levels.  This  was  done  to  facilitate  selection  of  a 
suitable  mesh  spacing  for  the  optimisation  algorithm  so  that  the  entire  crack-related  heat- 
production  zone  could  be  encompassed  within  the  10x10  computational  grid  used  in  the 
analysis. 

Fortunately,  thermography  affords  a  number  of  means  by  which  to  map  the  extent  of 
cyclic  plastic  deformation  in  a  material.  An  obvious  approach  is  to  apply  the  optimisation 
scheme  to  a  growing  region  of  data  centred  on  the  crack  tip  with  the  expectation  that  an 


7 


DSTO-TR-1190 


Figure  1:  Schematic  of  experimental  set-up  used  for  validation  of  the  energy  computation 
scheme.  Note  that  the  Aluminium  rod  serves  as  a  conduit  for  the  flow  of  heat  from  a  large 
thermal  mass  (not  shown)  to  the  Aluminium  plate. 


Figure  2:  Contour  of  integration  for  the  evaluation  of  the  power  density  within  region  A 
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nominal  input  power  (a.u.) 

Figure  3:  Comparison  between  the  flux-based  and  inverse  estimates  for  the  heat  source 


Figure  Schematic  showing  type-A  (left)  and  type-B  specimen  geometries. 


expansion  of  this  region  beyond  that  which  fully  contains  the  heat-production  zone  would 
produce  no  further  increase  in  estimated  total  power.  Although  likely  to  be  successful, 
the  approach  is  time  consuming.  A  more  expedient  alternative  is  to  presume  that  heat 
production  occurs  wherever  a  significant  temperature  increase  is  recorded.  This  of  course 
ignores  the  effect  of  heat  conduction,  which  will  tend  to  overstate  the  actual  zone-size. 

The  method  adopted  here  was  to  use  the  second  harmonic  thermal  response,  which 
is  that  component  of  the  measured  thermal  response  that  occurs  at  twice  the  loading 
frequency.  The  reason  for  its  use  is  two-fold.  The  process  of  mapping  the  second  harmonic 
response  is  simple  and  relatively  fast,  usually  achieved  in  under  a  minute.  And  secondly, 
the  technique  is  likely  to  be  effective  since,  due  to  its  inherent  non-linearity,  cyclic  plastic 
deformation  should  manifest  a  much  higher  second  harmonic  response  (as  well  as  higher 
harmonics)  than  regions  undergoing  purely  elastic  loading. 

Measurement  of  the  second  harmonic  component  was  achieved  by  performing  at  each 
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Figure  5:  Variation  in  measured  second  harmonic  response  due  to  the  application  of  an 
overload.  Crack  propagation  direction  is  from  left  to  right. 


pixel  location  a  real-time  cross-correlation  between  the  measured  infrared  response  and  the 
square  of  the  load  signal.  At  a  12.5  Hz  loading  frequency,  the  second  harmonic  response 
occurs  at  25  Hz  which  coincides  with  the  camera  sampling  frequency  and  would  result  in 
severe  aliasing.  To  avoid  this,  loading  was  applied  at  a  lower  frequency  of  10  Hz.  Since 
the  phase  difference  between  the  temperature  and  strain  evolutions  is  a  priori  unknown, 
both  in-phase  and  qnadrature  components  of  the  correlation  were  obtained,  yielding  a 
resultant. 


=  ^/Oii  +  olj)  (30) 

where  I  and  O  are  respectively  the  in-phase  and  quadrature  components  of  the  cross¬ 
correlation.  Figure  5  shows  a  sequence  of  maps  of  the  resultant  acquired  in  relation  to  an 
experiment  where  a  sample  exposed  to  persistent  cyclic  loading  was  suddenly  subjected 
to  an  overload. 
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The  three  images  {Figure  5)  relate,  respectively,  to  a  time  (a)  prior  to  application  of 
the  overload,  (b)  immediately  following  the  overload,  and  (c)  at  a  stage  where  the  crack 
had  grown  a  further  2  mm  in  length.  Variations  in  the  size  of  the  cyclic  plastic  zone,  as 
inferred  from  the  second  harmonic  response  measurements,  are  broadly  consistent  with 
the  effects  understood  to  follow  the  application  of  overload.  The  reduction  recorded  from 
images  (a)  to  (b)  for  example  are  consistent  with  one  or  a  combination  of  the  following 
overload  induced  mechanisms  (see  [11]):  compressive  residual  stress,  crack  tip  blunting 
and  crack  closure.  The  subsequent  expansion  in  zone  size  from  (b)  to  (c)  is  consistent 
with  the  crack  having  grown  through  the  overload  affected  zone. 

Based  on  these  second  harmonic  response  measurements,  an  analysis  window  of  11  mm 
X  11  mm  was  deemed  sufficiently  large  to  encompass  the  anticipated  cyclic  plastic  zone 
size,  and  accordingly  was  used  for  all  subsequent  experimental  work. 


4  RESULTS 

4.1  Specimen  type- A 

Seven  single-edge  notched  specimens  were  fatigue-tested  under  the  loading  regime  de¬ 
scribed  earher.  Temperature  response  data  were  acquired  at  irregular  but  frequent  inter¬ 
vals  during  the  fatigue  test  and  the  energy  dissipation  calculations  performed  off-line  at 
the  completion  of  the  testing  program.  Figure  6  shows  the  evolution  in  computed  power 
as  a  function  of  crack  length  for  all  of  the  tested  samples.  The  upward  monotonic  trend 
evident  here  is  intuitive  and  consistent  with  the  intensification  and  enlargement  of  the 
region  of  cycHc  pleistic  deformation  which  is  expected  as  the  stress  amplitude  increases 
with  crack  growth  imder  constant  amplitude  loading. 


Figure  6:  Dissipated  power  as  a  function  of  crack  length  for  type- A  geometry 


More  noteworthy  is  the  apparent  linear  relationship  between  the  power  and  crack 
growth  rate  {Figure  7).  Several  earher  studies  that  have  explored  the  relationship  be¬ 
tween  energy  and  fatigue  in  metals  (see  for  example  [9])  have  raised  the  possibility  that 
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structural  failure  of  an  element  of  material  may  occur  when  the  stored  energy  reaches  a 
certain  ‘critical’  level  characteristic  of  the  material.  If,  in  the  present  study,  cyclic  plastic 
deformation  occurred  with  a  constant  dissipation  function,  then  the  observed  linear  trend 
would  support  this  assertion.  Since  there  was  no  measure  of  the  applied  mechanical  work 
here,  no  assessment  of  the  dissipation  function  could  be  made,  and  so  the  link  is  largely 
speculative.  It  does  however  encourage  further  investigation. 

Considering  once  again  the  result  shown  in  Figure  7  it  is  notable  that  an  extrapola¬ 
tion  of  the  best-fit  trend  passes  almost  precisely  through  the  origin.  Prom  an  intuitive 
standpoint  this  seems  reasonable  in  suggesting  that  crack  growth  cannot  occur  without 
energy  dissipation.  Again,  it  has  not  been  possible  to  experimentally  verify  whether  the 
relationship  does  indeed  hold  as  the  crack  growth  rate  approaches  zero.  Attempts  were 
made  to  do  this,  however  it  became  apparent  that  a  much  improved  detector  sensitivity 
would  be  required  for  success  in  testing  imder  these  conditions. 


Figure  7:  Dissipated  power  as  a  function  of  crack  growth  rate  for  the  type-A  specimen. 
Error  band  corresponds  to  one  standard  deviation  of  the  data, 

It  is  worthwhile  commenting  on  the  increasing  scatter  as  a  function  of  crack  growth 
rate.  This  is  presumed  to  arise  in  part  as  a  result  of  the  increased  uncertainty  in  the 
measurement  of  crack  length  with  increased  rate  of  crack  growth. 


4.2  Specimen  type-B 

The  type-B  specimen  design  was  conceived  to  provide  a  simple  means  to  test  for  the 
effect  of  specimen  geometry  on  the  previously  established  linear  relationship.  The  results 
{Figure  8)  indicate  a  marked  deviation  in  the  evolution  of  the  energy  dissipation  rate  as  a 
function  of  crack  length  over  that  recorded  for  the  type-A  configuration.  A  reduction  in  the 
net  cross-sectional  area  is  expected  to  lead  to  higher  crack-tip  stresses,  and  consequently 
to  greater  crack-tip  plasticity  at  most  crack  lengths,  which  is  supported  by  the  data. 
Despite  the  clear  distinction  here  between  the  two  sets  of  results,  the  variation  in  the  rate 
of  heat-dissipation  with  crack  growth  rate  {Figure  9)  appears  largely  unchanged. 
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Figure  8:  Dissipated  power  as  a  function  of  crack  length  for  type- A 


and  type-B  specimens. 


Figure  9:  Dissipated  power  as  a  function  of  crack  growth  rate  for  type- A  and  type-B 
specimens. 
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4.3  Specimen  type-A  with  Overload 

Load  interaction  effects  caused  by  complex  loading  axe  known  to  play  an  important 
role  in  dictating  crack  growth  behaviour.  In  a  preliminary  attempt  to  assess  the  impact 
of  such  events  on  the  relationship  between  heat  dissipation  and  crack  growth  rate,  tests 
were  conducted  on  the  type-A  geometry  where  a  constant  amplitude  loading  sequence 
(described  earher)  was  augmented  with  a  single  overload  event  applied  at  a  nominal  crack 
length  of  15  mm.  Prior  experimentation  had  revealed  that  a  static  load  30%  higher  than 
the  peak  cyclic  load  would  lead  to  the  formation  of  a  substantial  monotonic  plastic  zone 
surrounding  the  crack  tip  (see  Figure  10),  and  consequently  facilitate  crack  closure  and 
the  formation  of  a  zone  of  compressive  residual  stress. 


Figure  10:  Monotonic  plastic  zone  size  produced  by  a  30%  overload  event. 

Figure  11  demonstrates  an  immediate  and  significant  impact  relating  to  the  overload 
event.  Initially,  the  energy  dissipation  rate  is  suppressed  relative  to  that  recorded  under 
benchmark  conditions,  maintaining  a  steady  value  as  a  function  of  crack  length  for  approx¬ 
imately  3  mm  of  crack  growth.  This  is  followed  by  a  period  of  rapid  ascent  in  dissipation 
rate,  which,  by  the  end  of  testing,  slightly  exceeds  that  recorded  in  the  absence  of  over¬ 
load.  Figure  12  compares  the  dissipation  rate  to  the  rate  of  crack  growth.  The  increase  in 
scatter  associated  with  the  computation  of  crack  growth  rate,  combined  with  the  relatively 
small  data  set  limits  the  scope  of  comment  possible  here.  Although  the  data  does  suggest 
that  load  interaction  effects  are  likely  to  play  an  important  role  in  this  relationship,  it  is 
difficult  to  speculate  on  the  physical  mechanisms  at  work.  For  instance,  the  presence  of 
compressive  residual  stress  is  known  to  supress  the  rate  of  crack  growth  to  levels  below 
that  predicted  by  a  stress-intensity  factor  calculation,  leading  to  so-called  anomolous  be¬ 
haviour.  Measurements  of  energy  dissipation  however  should,  unlike  the  stress  intensity 
factor,  reflect  the  presence  of  a  compressive  residual  stress  since  the  level  of  cyclic  plastic 
deformation  is  also  suppressed,  in  essence  leading  to  some  level  of  cancellation.  Clearly, 
this  is  mere  speculation  and  cannot  be  tested  with  the  available  data,  but  is  nevertheless 
an  interesting  possibility  worth  further  investigation. 
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Figure  11:  Dissipated  power  as  a  function  of  crack  length  for  the  type- A  configuration  with 
and  without  overload. 


Figure  12:  Dissipated  power  as  a  function  of  crack  growth  rate  for  the  type- A  configuration 
with  and  without  overload. 
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5  CONCLUSION 

A  numerical  scheme  has  been  presented  to  allow  for  the  remote  assessment  of  the  energy 
dissipation  associated  with  crack-growth  in  metals  on  the  basis  of  remote  radiometric 
temperature  measurement.  After  providing  experimental  verification  of  the  scheme,  it 
was  applied  in  an  investigation  of  the  relationship  between  energy  dissipation  and  the  rate 
of  crack  growth  in  mild  steel  exposed  initially  to  consteint  amplitude  cyclic  loading.  A 
hnear  relationship  was  foimd  and  then  reaffirmed  in  tests  on  a  more  complex  geometry, 
indicating  some  level  of  geometry  independence,  an  encouraging  result  from  the  standpoint 
of  developing  a  useful  crack  characterisation  parameter.  The  relationship  did  not  appear 
to  be  maintained  when  load  interaction  effects  were  included  by  means  of  a  static  overload, 
suggesting  that  other  factors  not  considered  in  the  present  study  are  involved. 
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